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Abstract. In this paper we illustrate how non-stochastic (max,+) techniques 
can be used to describe partial synchronization in a Discrete Event Dynamical 
System. Our work uses results from the spectral theory of dioids and analyses 
(max,+) equations describing various synchronization rules in a simple network. 
The network in question is a transport network consisting of two routes joined at 
a single point, and our Discrete Events are the departure times of transport units 
along these routes. We calculate the maximum frequency of circulation of these 
units as a function of the synchronization parameter. These functions allow us 
further to determine the waiting times on various routes, and here we find critical 
parameters (dependent on the fixed travel times on each route) which dictate 
the overall behavoiur. We give explicit equations for these parameters and state 
the rules which enable optimal performance in the network (corresponding to 
minimum waiting time). 



1 Introduction 

Considerable advances have been made in recent years in the algebraic descrip- 
tion of Discrete Event Dynamical Systems (DEDS). Chief among these is the 
approach using the so-called (max,+) semiring [Q, with corresponding develop- 
ment of algorithms for spectral calculations M and concrete applications H . 

In this paper we illustrate how non-stochastic (max,+) techniques can be 
used to describe partial synchronization (in a manner to be defined below) in 
a DEDS. Our work uses results from the spectral theory of dioids and analyses 
(max,+) equations describing various synchronization rules in a simple network. 

The paper is structured as follows. In Section || we review necessary math- 
ematical properties of the (max,+) semiring. Section ^ introduces the network 
on which our DEDS is based, and in Section [| we describe the different (max,+) 
equations we will solve. The resulting solutions lead, in Section ^|, to critical 
phenomena in the network, and we calculate which among the (max,+) equa- 
tions give optimal results. A numerical example illustrates our results in Section 
we conclude in Section fH. 



2 (max,+) semiring 



The (max,+) semiring (or dioid) 7^. max is the set 1Z U {—00} with the two op- 
erations ffi, © defined by a © b = max(a, 6) and a © & = a + b. Note that both 
operations possess identity elements, —00 for © and for ©, but while © is 
invertible, © is idempotent: a © a = a. In the remainder of this paper we 
will often omit the explicit multiplication symbols x,©, as is the convention. 
Our convention will be that if an expression contains any explicit symbol from 
we will assume all such hidden symbols are also from lZ max : a + bc means 
a+ (b x c) , a® be means a © (b © c) , and (a © b) (c + d) means (a © b) © (c + d) . 

These operations can be extended to matrices A,Be R™ by defining 
(A © B)ij = Aij © Bij and (A © B)ij = (Bk(Aik © -Bfcj)- As this may seem 
unusual to the novice reader, a small example is in order. 

/ (3©4)ffi(-l©0) \ / 7 
\ (0 © 4) © (5 © 0) / V 5 

We now look at solutions to the (max,+) spectral problem A © X = A © X. 
Theorem Let A e • There exists a maximal eigenvalue 

n 

A = (A tll2 ©A l2l3 ©...©A (fcll )/fc. (2) 

If we extend our notation to a n d = a©a© . . .©a = n x a and Tr(A) = ©fc = i ^4fefe 
equation (j|) reads 

n 

A = 0(Tr( J 4 fe )) 1 / fc . (3) 
fc=i 

In the standard way, we associate with A a weighted digraph (precedence graph) 
Q{A) with nodes M = {1,2, ,..,n] and edges £ = {(i,j)|A,-j 7^ -00}. The 
weight of edge is simply Aji. Denoting by i — > j the existence of a path 
from i to j we define A to be irreducible (and Q(A) to be strongly connected) 
iff i — ► j V 1 < i, j < n. 

Theorem 2.Q If A E R™ is irreducible, then X as given by (g) is unique. 

Note that from a graph-theoretic point of view, A is merely the maximum cycle 
mean of G(A). 

Our DEDS will be described by equations in TZ m &x'- Specifically, denoting 
the time of the fcth occurence of event i by Xi(k), we are interested in equations 
of the form 




2 



(4) 



where we have a finite family of 1 < i, j < n events and 1 < I < k — 1. For 
linear equations in 7£ max , (Q) becomes 

a*(fc)=©(/ii®ay(*-0) (5) 
i,* 

and we will show in Section |] how this can be written as a matrix equation with 
unitary retard 

Xi(k) = 0(A y - ® ^-(fc - 1)). (6) 
j 

We are interested in the asymptotic properties of this system as k — > oo so we 
study limfc^oo A k . 

Theorem 3.Q0][^| If A G ^™iax * s irreducible, there exists integers M,c(A) 
such that 

A k®c(A) = ( X y(A) ^ A k y fc > M. (7) 

We call c(A) the cyclicity of A (sec |l[ for details of how to calculate this 
integer). 

3 Model 

We consider an elementary model of two tourbuses with partial synchronization. 
Each tourbus goes on a fixed circular route R\ and i?2 , stopping every so often 
to drop off and pick up passengers. The tour buses meet at one location, a 
downtown station S, where passengers can pass from one bus to the other one. 
The physical network consists of two loops meeting at S, as shown in Figure 1. 

We assume fixed travel times T\ and Ti for R\ and R2 respectively, and 
without loss of generality let T\ < T2. The simplest model is to describe this 
system using the (max,+) equations X(k + 1) = A ® X(k), with 

A=( ai1 ai *),X(k)=( ) (8) 

V fl 21 022/ V X2 \ k > J 

and Xi{k) is the departure time of bus i from station S along route We 
consider the following extreme cases. 
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Figure 1: Physical network of routes Ri,B,2 linked at station S. 



CASE Ml (no synchronization): 

We have 




oo T 2 



— oo 



) 



(9) 



so that the graph Q{A) is disconnected into two parts each having its own cycle 
mean Ai = T\ and A2 = T 2 as shown in Figure 2 (the nodes here correspond to 
departure times of buses). 



Figure 2: Precedence Graph. One bus each on R\ and R 2 . Not Synchronized. 

CASE M2 (complete synchronization): 

We have 



so that the graph Q{A) is strongly connected. There are 3 circuits in this graph 
with cycle means T 1; T 2 , and Tl ~ 2 ~ T2 , so the maximum cycle mean gives an 
eigenvalue of A = T 2 (see Figure 3). 

4 Partial Synchronicity 

For the case of complete synchronization in section ||, it is clear that the large 
value of A = T 2 slows down the overall system. We further suppose in this 




T 2 




(10) 
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Figure 3: Precedence Graph. One bus each on R x and R2. Completely Syn- 
chronized. 



section that we cannot alter the travel time T2, but that we have a large amount 
of additional stock (buses) which we can add to the network. 

We add to route R2 a large number m of buses such that the interval between 
departure times on this route becomes a very small number S (which we can 
think of as T2/1T1 if required). Denote by y\(k), (k), ... , y m {k) the departure 
times of the m new buses, which operate according to the rules 

Vi(k)= 5<8>2/i_i(fc-l) 

Vo(k)= 8®y m (k-l). (11) 

For this route we do not produce a timetable, as we assume to a first approxi- 
mation there is virtually no waiting for a bus. 

For this system, it would seem there should be no interdependence between 
X\(k) and X2{k). However, what of a passenger who wants to make a non-stop 
round trip R\ — > R2 — > R\ - would some synchronization rule benefit such a 
passenger? To investigate further we set T2 = nT x + r and look at the 7?. max 
equation x x (k) = f(x x (k — I)) for different values of I and linear /. 

CASE PI: 

Xl {k) = T lXl (k - 1) © T x T 2 x x {k -(1 + 1)) (12) 

Using n artificial states (i.e. states that do not in fact correspond to the depar- 
ture of any physical bus) zi, Z2, ■ ■ ■ , Zi, we can rewrite ([12]) as 

xi(k) = Tixx(k- 1) ®rT x zi{k- 1) 

Zi(fc) = T? /l ®z i - 1 (k-i),i = l,...,l (13) 

Zo(fc) d = xi(k). 

At the expense of increasing the dimension, our system of equations now has 
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unitary retard as in (^|) with column vector 



X(k) = 



( Mb) \ 

zi{k) 
zi-i(k) 



(14) 



and matrix 



.4 



— oo 
— oo 

V nTxjl 



-oo 



-oo nTi/l 



-oo nTi/l 



-oo — oo — oo 



— oo \ 

— oo 
— oo 

— oo J 



The full precedence graph Q{A) for the system x\, z%, . . . , Zi, X2, j/i, ■ ■ ■ , y m is 
shown in Figure 4. The disconnected parts of the graph G(A) have maximum 




Figure 4: Precedence Graph for CASE PI, equation (|12| 
cycle means A2 = S and 

r Ti if Z > n 

Al = < (n+l)Ti+r- 

L z+i 



(15) 



otherwise 

since Ti©r = T% (by definition of r) and, for Z > n (equivalently Z = Z©(n+1)), 



(n+l)Ti+r (ra+l)Ti+Ti _ (n+2)Ti <^ rp 
<v , , -, — J+1 ^ -11. 



J+1 



(16) 



Note that Z > n effectively corresponds to no synchronization, while for the 
choice Z < n, (n — I + 1) out of every (n + 1) buses are synchronized. The 
minimum possible Ai is Ai = T\ given by Z > n. 

CASE P2: 



x 1 {k)=T 1 x 1 {k-l)®rT[ +l x 1 {k-{l + l)) 



(17) 
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As in the previous case PI, we can rewrite (O) as 



x x {k) = T 1 x 1 {k-l)®rT 1 zi{k-l) 
Zi(k) = T 1 ®Zi-i{k-l),i = l,...,l 
z {k) = xi(k). 



(18) 



and the system is described by X(k + 1) = A ® X(k) with X(k) as in (fTl) and 



A = 



I T\ T\ + r — oo — oo 
— oo — oo T\ — oo 
— oo — oo — oo T\ 



-oo \ 
-oo 

- 30 



y Ti — oo — oo — oo • • • — oo J 
From the precedence graph shown in Figure 5 we have A2 = 5 and 

(l + l)Ti+r 



Ai = 



1 + 1 



since obviously 



(I + l)Ti + r 



(19) 



(20) 



Z + l 

for all values of /. This system corresponds to one out of every (I + 1) buses 




Figure 5: Precedence Graph for CASE P2, equation ( |l~7|) 

being synchronized. The limiting case is 

lim (/ + l)T 1 + r =ri 

i^oo / + 1 



(21) 



so the minimum value of Ai occurs at I = 00. 
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5 Optimal Solutions 



In both cases PI and P2 of Section || the minimum value of Ai is T\. In PI it 
is given by any choice I > n while in P2 it is given by the extreme value I = oo. 
These values both correspond to no synchronization, and this is the result we 
would expect, that the system runs faster if the buses do not wait for each other. 

If our goal is to choose the synchronization rule which maximizes the speed 
of the system, the problem is solved. In this Section we look at the separate 
problem of minimizing the waiting time of passengers that make a circular 
journey over all or part of the network (by circular we mean a journey that 
begins and ends at the same point). We itemize such journeys Ji as follows: 



Jl 


Ri - 


+ R-y 


Ji 


Ri - 




h 


Ri - 


*Rx 


Ja 


Ri - 


*Rl 



Ri 
R2 



For cases PI, P2 we look at regions I < n and I > n to give synchronization 
rules Si as follows: 



Equation (12), I < n 
Equation Jl2]), I > n 
Equation I < n 
Equation dlTl) ,1 > n 



St 
S 2 
S 3 
S4 

For i,j e {1,2,3,4} we pair each journey Ji with rules Sj to give 16 differ- 
ent models. Let Wij be the total waiting time on journey Ji with rule Sj. 
These waiting times can be broken into two components, the waiting time be- 
fore boarding the bus w\j and the waiting time in transit w\p which can only 
occur for journeys Ji and J3. From our calculation of Ai and A2 in Section || 
we determine w\j , w\j and = w\a + wjj to be the values given by equations 
@-@ in the Appendix. 

We now make some observations on the solutions given by equation (|3l|). 
Our goal is to minimize the waiting times, or more specifically the average 
waiting time per passenger as a function of the synchronization parameter I: If 
n, people take journey Ji and N = Uj then we want to find the minimum of 



IF, 



dcf J2j n j w ji 



N 



(22) 



Note that 



• wn is minimized by choosing the largest possible value of I, i.e. I = n. 

• w i2 = lim^oo w i4: 

• similarly is minimized by choosing / = n. 
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Finally examine wa. Wu and W34 are both clearly minimized by I = 00. To 
minimize W24 we view it as a continuous function of Z and solve 



dw 2 4 d (3Z - 2n + 1)7\ + (2n - 2/ + 3)r 
dZ ~ <tt 2(i + 1) ( J 

Z drops out of equation (p3|), meaning the existence of an extremum does not 
depend on I but rather on the solution of (2n + 2)7\ = (2n + 5)r. Thus we have 
a critical point 



/2n 



V 2n 



Tl (24) 



The slope of the function W24,(l) can now be written in terms of the critical 
remainder r c as 

(hv24_ _ (2n + 5)(r e - r) 
dZ " (Z + l)2 • 

Hence W24Q) is 

• strictly increasing for r < r c : minimum is obtained by choosing lowest 
value of Z, i.e. I = n + 1 

• strictly decreasing for r > r c : minimum is obtained by choosing highest 
value of Z, i.e. Z = oo 

• constant for r — r c : choose any value of Z 

In equation ([}2|) we rewrite W with the optimal values of Z chosen as indicated 
above: for we insert in the matrix W n (= W min ) the values at Z = n + 1, 
since the values at Z = 00 correspond to column 2, Wi2- 

Now that we have chosen Z we should choose the synchronization rule Sj that 
minimizes Wi. In rows 1, 3 and 4 of W m it is clear that the minimum values 
are in column 2 corresponding to 52. Let us look at row 2, w^J. In Figure 6 we 
plot wl 7 !, an d W24 as functions of r. 



2n 



T x (26) 



2n + 3, 

along with equation ([24|) defines three different regimes as follows: 

• r > r s : W22 < W21 < ^24 

• r c < r < r s : w 2 \ < wV^ < 

• r < r c : w 2 ™ < < 

Further we note the importance of the critical value r s is that for r > r s rule 
5*2 gives shortest waiting times for all journeys Ji, and in this regime W2 is 
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TV2+Ti/(n+2)~ 
Ti a 



Figure 6: w^l,w^ and w™ 4 plotted against < r < T\ (see equation (|32|)): 
The critical points are r c = (2n + 2)T 1 /(2n + 5) and r s = (2n + 2)T 1 /(2n + 3). 



minimal irrespective of the values rij (see equation (|22|)). For r < r s the values 
of rii dictate which Wi is minimal: In particular, if 



n 2 > 



(2n + 3)(r s - r) 



(3n-i + n 3 ) 



(27) 



then Wi is minimal. We now have a rule that governs optimum performance in 
our network: If r > r s implement rule S2, otherwise monitor passenger numbers 
m and if at any stage equation ( p7j ) holds, implement rule Si. 

Observe further from equation ( p6| ) that 4Ti/5 < r s <T\, since 1 < n < 00. 
This means that even for a random choice of T\ and T2, the chances that r is 
less than r s are greater than 80%: 



P(r < r s \ arbitrary Ti,T 2 ) > 0.8 



(28) 



6 Example 



Let Ti = 3 units of time and n — 2 so that r s = 18/7 as given by equation (f26|). 

To examine firstly the regime r > r s , let r = 8/3, given by Ti = 26/3 
units of time. Under rule S2 we calculate w™ 2 = 33/18 with actual departure 
times x\{k) — 0,3,6,9,12,15,... Under rule S3 we have = 35/18 with 
departure times xi(k) = 0, 3, 6, 35/3, 44/3, 53/3, . . . Thus rule S2 gives optimal 
performance on J 2 - 

Secondly, let r — 7/3 given by T 2 = 25/3, so that we are in the regime 
r < r s . Rule ^2 gives = 39/18 with departure times Xi(k) = 0, 3, 6, 9, 12, 
15,. . . Under rule £3 we have = 34/18 with departure times x\{k) = 0, 3, 
6, 34/3, 43/3, 52/3,. . . Here optimal performance on J 2 is given by rule S3. 

Note in both regimes of this example the departure times x\{k) under S 2 
are identical, the reason for the differences in wW lying in the transfer waiting 
time. 
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7 Conclusion 



In this paper we have solved the (max,+) equations for a specific model with 
two routes joining at one point. We have shown how partial synchronicity can 
be described by these equations by introducing a synchronization parameter I. 
Using results from (max,+) spectral theory we have calculated the maximum 
frequency (minimum cycle mean) of buses on such routes. These numbers al- 
low us to calculate waiting times on various routes: We observe new critical 
behaviour which depends on the relation between T\ and T 2 , the (fixed) travel 
times on each route. Finally we state the rules for optimal performance in this 
network. 

Further work is envisaged in looking at more general (max,+) equations, in 
particular a system with m\ buses on R\ and ni2 on i?2- Another direction in 
which this work could possibly be extended is to have further routes R3, R4, . . . 
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Appendix 
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(29) 



(30) 
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